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ABSTRACT 



Context. The NASA exoplanet search mission Kepler is currently providing a wealth of light curves of ultra-high quality from space. 
Aims. We used high-quality Kepler photometry and spectroscopic data to investigate the Kepler target and binary candidate 
KIC 5988140. We aim to interpret the observed variations of KIC 5988140 considering three possible scenarios: binarity, co-existence 

■ of 6 Set- and y Dor-type oscillations, and rotational modulation caused by an asymmetric surface intensity distribution. 
Methods. We used the spectrum synthesis method to derive the fundamental parameters T cff , logg, [M/H], and vsin; from the newly 
obtained high-resolution, high S/N spectra. Frequency analyses of both the photometric and the spectroscopic data were performed. 

^ C"| Results. The star has a spectral type of A7.5 IV-III and a metallicity slightly lower than that of the Sun. Both Fourier analyses reveal 

O [ the same two dominant frequencies F 1 =2F2=0.688 and F2=0.344 d -1 . We also detected in the photometry the signal of nine more, 

significant frequencies located in the typical range of 6 Set pulsation. The light and radial velocity curves follow a similar, stable 
double- wave pattern which are not exactly in anti-phase but show a relative phase shift of about 0. 1 period between the moment of 
minimum velocity and that of maximum light. 

Conclusions. Such findings are incompatible with the star being a binary system. We next show that, for all possible (limit) con- 
figurations of a spotted surface, the predicted light-to-velocity amplitude ratio is almost two orders larger than the observed value, 
which pleads against rotational modulation. The same argument also invalidates the explanation in terms of pulsations of type y Dor 
(i.e. hybrid pulsations). We confirm the occurrence of various independent 6 Set-type pressure modes in the Kepler light curve. With 
respect to the low-frequency content, however, we argue that the physical cause of the remaining light and radial velocity variations 
1 of this late A-type star remains unexplained by any of the presently considered scenarios. 

Key words. Asteroseismology - Kepler - Stars: variable: <5 Set - Stars: atmospheres - Stars: abundances - Stars: rotation - Stars: 
individual: HD 188774 
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,— 1 . 1. Introduction tween 6 and 15 observed by Kepler has been performed by 

, xr . OA , , ,. „™„ lUvtterhoeven et al. I (1201 lbl) . In this study. 63% of the sample 

■ The NASA exoplanet search mission Kepler, launched in 2009, was asgi d tQ Qne of ^ main dasses of lsati A-F type 

L; ; is currently changing our views about pulsating stars and aster- stars . 2?% arg classified as § Sct stars (206 stars)> 23% as h brid 

£ ■ oseismology thanks to the collection of light curves on all kinds lsators exhibiting both 6 Sct- and y Dor-type oscillations (171 

^ ■ of variable stars. The nearly continuous time series with micro- g(arg)> and u% afj Dor ^ (m ^ The remaini stars 

»H . magnitude precision opens up opportunities for high-quality were classified as rotationa n y mo dulated/active stars, binaries, 

CO . and in-depth astero s eismic studies with unprecedented detai l or ^ ^ show nQ ^ iodk variab0it xhe majorit of 

jJenkins et aL 1 2010| ; | Chaphnetal. || 20lU | Beddingetal. || 20lU ±e hyhM pulsators shows f reque ncies with all kinds of period- 

| Beck et al. || 2011|) . icities within the y p o r and 5 Sct range, which is a challenge 

A first general characterization of the pulsational be- for the current models. |Ujlterfioe^ri^ay also found 

haviour of 750 candidate A-F type stars with magnitudes be- indications for the 

existence of 6 Sct and y Dor stars beyond the 



* Based on data gathered with NASA's Discovery mission Kepler and . , ,. , rrr-. — ■ — r nr^Tm. .I^r 2 — , , n 

with the Hermes spectrograph, installed at the Mercator Telescope, op- inthe studies by| Gngahcene et al. | (|2010|) and| Tkachenko et al. | 

erated on the island of La Palma by the Flemish Community, at the 



edges of the cu rrent theoretical instabilit y str ips, as confirmed 
in the si 
J20T2I) . 



Spanish Observatorio del Roque de los Muchachos of the Instituto de 
Astroflsica de Canarias, and with the 2-m Alfred- Jensch telescope of 

the Thiiringer Landessternwarte Tautenburg. Ground-based follow-up observations of the Kepler targets 

** Table A.l is available in electronic form only at the CDS have been organized with the goal to obtain precise fundamental 

via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via paramete rs needed for the seismic m odeling of pulsating stars 

http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-i-A/ (see e.g. lUvtterhoeven et al. I [20 1 lal) . A first sub-sample of A- 
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F type stars has been recently presented by ICatanzaro et al. I 
(1201 lh . 

We focus here on the Kepler target KIC 5988 140 
(HP 188774, Keple r magn itude of 8.852). In the study by 
lUvtterhoeven et al~1 (1201 lbl) . this late A/early F-type star is 
classified as a binary with 6 Set pulsations. At least 12% of the 
entire sample was identified as belonging to a binary or a multi- 
ple system. We further remark that its light curve resembles that 
of KIC 9664869 (see panel c of Fig. 5 in lUvtterhoeven et al. I 
1201 lbl) which is in turn classified as a hybrid pulsator of the mid- 
dle type, i.e. neither of dominating 5 Set type nor of dominating 
y Dor type. ICatanzaro et al. I d2011l) investigated KIC 5988 140 
based on low and medium resolution spectroscopic data and 
reported the fun damental atmospher ic param eters shown in 
Table Q] Unlike lUvtterhoeven et al. I d2011bl) . these authors 
concluded that KIC 5988140 belongs to the hybrid 6 Set - 
y Dor class. In this study, we use new high-resolution spectra 
of KIC 5988140 to investigate its radial velocity (RV) curve and 
to re-determine its fundamental parameters, with the goal to 
identify the possible cause(s) of the observed photometric and 
spectroscopic variations. 



2. Observations 

Our analysis is based on high-quality photometric and spectro- 
scopic data. The light curves have been gathered by the Kepler 
satellite in the so-called Long Cadence mode (with a time reso- 
lution of 29.4 min). The Kepler data are made available in quar- 
ters, i.e. periods between two spacecraft rolls. The (uncorrected) 
data of different quarters have different zero-points and suffer 
from instrumental drift which may vary from quarter to quar- 
ter. To correct for these, we used a package of Fortran routines 
(kindly made available by L. Balona), which removes the jumps 
between the quarters and the drifts in an automatic way. The pro- 
cessed data represent normalized relative intensities expressed 
in parts per thousand (ppt). Using this procedure may affect 
some of the low est frequencies found in a period-search analysis. 
lBalonald2011l) tested this assumption and concluded that the fre- 
quencies above 0.1 d~' are not affected, while the amplitudes of 
those below 0. 1 d~' may change slightly. For this case study, we 
analysed the Kepler data from seven different quarters, i.e. Q0- 
Q2,Q4-Q6,Q8, which resulted in 22234 measurements in total 
and from which we extracted 20668 useful measurements. The 
Nyquist frequency is 24.47 d _1 and the fr equency resolution is 
0.002 2 d -1 (for the time base of 682 days) dLoumos & Deeming I 
Il978l) . Two small aliasing peaks located at respectively 0.0033 
and 0.0052 d _1 , and caused by the missing quarters Q3 and Q7, 
can be identified in the corresponding spectral window. 

We also acquired a series of high-resolution spectra in two 
consecutive years. In 2010, the spectra were taken with the 
Hermes spectrograph attached to the 1 .2-m Mercator telescope 
(Observatorio del Roque de los Muchachos, La Palma, Canary 
Island). Hermes (High Efficiency and Resolution Mercator 
Echelle Spectrograph) is a fibre-fed spectrograph which sam- 
ples the entire op tical range from 38 to 900 nm with a reso- 
lution of 85 000 dRaskin et al. 11201 lb . The spectra acquired in 
2011 have been taken with the Coude-Echelle spectrograph at- 
tached to the 2-m telescope at the Thiiringer Landessternwarte 
Tautenburg. The Tautenburg spectra have a spectral resolution 
of 32 000 and cover the wavelength range from 472 to 740 nm. 
Exposure times of the same order as the temporal resolution of 
the Kepler data were used in both cases (i.e. between 6 and 30 
min depending on the conditions). 



Table 1. Fundamental parameters and the 1-cr formal errors as derived 
from the mean spectrum. The KIC parameters and their typical errors 
(in parentheses) are listed for comparison. 



[M/H] 


T eff (/0 


logs 


v sin £ 


SpT 


(kms ') (kms ') 



this work 

-0.30(0.05) 7600 (30) 3.39(0.12) 52.0(1.5) 3.16(0.20) A7.5IV-III 

KIC 

-0.54(0.50) 7451 (200) 3.54(0.50) — A8IV-III 



Catanzaro et al. (2011) 
7400 (150) 3.7 (0.3) 



The data have been reduced using the dedicated pipeline 
for Hermes spectra and standard ESO-MIDAS packages for the 
extracted, merged Hermes spectra and the Tautenburg spectra, 
respectively. The data reduction included bias and stray-light 
subtraction, cosmic rays filtering, flat fielding, wavelength cal- 
ibration by ThArNe lamps, order merging, and normalization 
to the local continuum. All the spectra were additionally cor- 
rected in wavelength for individual instrumental shifts by using 
a large number of telluric O2 lines. The cross-correlation tech- 
nique was used to estimate radial velocities (RVs) of the individ- 
ual spectra. All the spectra have been shifted in wavelength ac- 
cording to their derived RVs and co-added to build two high(er) 
signal-to-noise ratio (S/N) mean spectra. The Hermes (respec- 
tively Tautenburg) mean spectrum was next analysed individu- 
ally. Since both analyses gave consistent results, and because 
the Hermes mean spectrum covers a wider spectral range and is 
of excellent quality (S/N = 180-200 at A ~ 500 nm), we will re- 
fer to the latter (only) as 'the mean spectrum' in the following 
sections. 

The journal of observations is given in Table IA.1I The as- 
sociated errors quoted in Table IaTI are the 1-cr formal errors 
obtained during the fitting process. 

3. Spectroscopy 

3.1. Spectrum analysis 

The mean spectrum of KIC 598 8140 was analysed by m eans of 
the GSSP programme package dTkachenko et al. 1120121) which 
is based on the method of synthetic spectra. The idea is to use 
a large, dense enough grid of pre-computed theoretical spec- 
tra and perform a spectrum-by-spectrum comparison with the 
observations until a global minimum of x 2 m parameter space 
is found. In the first step, synthetic spectra are computed on 
a grid of effective temperature T e ff, surface gravity logg, pro- 
jected rotational velocity v sin i, microturbulent velocity and 
metallicity [M/H]. The derived metallicity is used as an initial 
guess for the chemical composition of the star and individual 
elemental abundances are iterated together with the other four 
fundamental parameters in the second step. The measurement 
uncertainties of T e ff, logg, vsin i, [M/H], and the individual 
abundances are determined as 1-cr confidence levels obtained 
from ^-statistics. For the calculation of synthetic spectra, we 
used the code SynthV dTsvmbal I IT996) based on the local ther- 
mal equilibrium (LTE) theory whereas atmosphere models have 
been compute d with the most rec ent, parallelized version of 
the LLmodels dShulvak et al. 1120041) . All synthetic spectra were 
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Table 2. Abundances for individual elements sorted according to their 
uncertainties. N is the number of lines with a line depth larger than rms. 
The uncertainties quoted between parentheses are given in terms of last 



digits. 














ilized flux 


PlpTTlPTlt 
J ^ 1 1 H_ J n 


N 


N I 


N II 


solar 


1A/H1 


dA 


Fe 


861 


753 


108 


^t.59 


^.99(05) 


-0.40(05) 


S 


Ca 


55 


42 


13 


-5.73 


-5.88(05) 


-0.15(05) 


5 


Ti 


162 


23 


139 


-7.14 


-7.42(09) 


-0.28(09) 




Cr 


149 


71 


78 


-6.40 


-6.60(11) 


-0.20(11) 




Mg 


32 


25 


7 


^t.51 


^.89(15) 


-0.38(15) 




Ni 


80 


74 


6 


-5.81 


-6.05(16) 


-0.24(16) 




Sc 


24 





24 


-8.99 


-9.17(19) 


-0.18(19) 




Mn 


42 


35 


7 


-6.65 


-7.07(21) 


-0.42(21) 




V 


44 


2 


42 


-8.04 


-8.27(23) 


-0.23(23) 




c 


14 


14 





-3.65 


-3.83(24) 


-0.18(24) 




Y 


28 





28 


-9.83 


-10.01(25) 


-0.18(25) 


X 


Si 


13 


4 


9 


-4.53 


^1.77(31) 


-0.24(31) 


3 

5= 


Co 


14 


13 


1 


-7.12 


-7.43(36) 


-0.31(36) 


u 


Sr 


4 





4 


-9.12 


-9.02(42) 


+0.10(42) 


N 


Ba 


7 





7 


-9.87 


-9.57(46) 


+0.30(46) 




Eu 


6 





6 


-11.52 


-11.47(49) 


+0.05(49) 


O 

c 



Note: Columns 1 to 4 represent the element designation, the total num- 
ber of spectral lines in the considered wavelength region, and in the first 
and the second ionization stage, respectively, while columns 5 to 7 list 
the solar composition, the individual abundances, and the deviation of 
those with respect to the solar values jAsplund et al. f 2009). 



computed in the 3800-5700 A wavelength region which is al- 
most free of telluric contributions but includes several hydrogen 
lines of the Balmer series and a large number of metal lines. A 
more detailed description of the method and its application to 
the spectra of Kepler ySCep and SPB candidate stars as well as 
6 Set and y Dor candidate stars are g iven by iLehmann et al. I 
d2011l) and bv lTkachenko et al. I (120121) . respectively. 

Table Q] gives the results of the spectrum analysis in com- 
parison with the Kepler Input Catalog (KIC) values. The errors 
of the KIC parameters are typical values, i.e. +200 K for the 
effective temperature and +0.5 dex for both the surface gravity 
and the metallicity. The derived fundamental parameters agree 
within the (KIC) errors, though there is a rather large difference 
in the value of [M/H] between ours and the KIC one. The new 
fundamen tal parameters also agre e with respect to the values re- 
ported bv ICatanzaro et alTI d201 lh . Spectral types and the lumi- 
nosity classes listed in Ta ble Q] were derived usin g interpolation 
in the tables published bv lSchmidt-Kalerl dl982l) . 

Table [2] lists the abundances of individual elements as de- 
termined from the mean spectrum. Apart from Sr, Ba, and Eu, 
the abundances of all elements are compatible with the derived 
metallicity of the star. The three elements mentioned above 
show considerable overabundances but their uncertainties are 
also large. 

Figure Q] shows parts of the observed spectrum (black line) 
in comparison with the two synthetic spectra computed respec- 
tively based on the parameters derived by us (dashed line) and 
on those listed in the KIC (dotted line). Obviously, the newly de- 
rived fundamental parameters provide a better fit of the observed 
spectrum in both illustrated wavelength regions. The same is true 
for the majority of the other metal lines found in the considered 
wavelength range. 
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Fig. 1. Fit of the observed spectrum of KIC 5988 140 (black) in two dif- 
ferent wavelength regions assuming T e ff = 7600/7400 K, log g = 3.4/3.5 
and [M/H] = -0.3/-0.5 dex (resp. dashed red and dotted green lines). A 
colour plot is provided in the online version. 



j Spectral window 




nJ 

DC ? 



Frequency (1/d) 

Fig. 2. Successive periodograms of the ground-based RV data of 
KIC 5988 140 (top: spectral window, bottom: after prewhitening for Fi 
and F 2 ). 



3.2. Radial velocity variations 

The RVs measured (by cross-correlation, see Sect. |2]i from 
all recently observed spectra (i.e. 40 measurements in to- 
tal) were searched for possible period i cities by means of the 
Period04 programme riLenz&Breger 1 120051) . Figure [2] illus- 
trates the Fourier spectrum of the RV data set. The corre- 
sponding spectral window reveals the presence of the 1 year 1 
as well as 80% high 1 day -1 alias peaks. The dominant pe- 
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Table 3. Frequencies detected in the RV data of KIC 5988 140 



..1 c o . 



. iz i r ' cnoo a An, 



ID 


Frequency 


Amplit. 


Phase 


S/N Note 




(± error) 


(± 0.05) 








d- 1 


km/s 


2;rrad 




Fi 


0.6880 (0.0038) 


2.93 


0.57 


33 = 2F 2 


F 2 


0.3440 (0.0002) 


1.34 


0.91 


15 


F 


0.0013 (fixed) 


1.41 


0.04 


15 



riod Pi is detected at 1.4534+0.0085 d. After prewhitening 
with the corresponding frequency ,Fi=0.6880±0.0038 d , a 
Fourier analysis of the residuals reveals a second significant 
frequency F2=0. 3440+0. 0002 d -1 corresponding to the period 
P 2 =2.9070+6.0002 d (see Table [3]). We computed 100 Monte 
Carlo simulations to estimate the errors on the frequencies and 
amplitudes (cf. Period04). This second frequency appears quite 
significantly after prewhitening, and is not affected by the alias- 
ing peaks located at ~0.5 d -1 in the spectral window (as will be 
confirmed in Sect.|4j>. We also remark that the ratio between both 
frequencies is exactly 2. The third frequency, F* = 0.0013 d _1 , is 
apparently caused by a small shift in zero-point between the RVs 
originating from two different instruments. We kept its value 
fixed during the Monte Carlo error computations. The tabulated 
ampli tude-to-noise ratios much larger than 4.0 dBreger et al. I 
1 19931) show that all three frequencies are significant. The frac- 
tion of the variance removed after prewhitening with three fre- 
quencies, (1 - R), equals 87%. 

Figure[3]illustrates the RV curve based on the Tautenburg and 
Hermes spectra, phased on the period P2 (top), as well as a small 
portion of the Kepler light curve (bottom). Phase zero in the RV 
curve corresponds to the time of the highest maximum observed 
in this part of the light curve which occurs at (B)JD 2455087.62. 
The following ephemeris has been adopted: 



(B)JD, 



maxJight/RV 



(B)JD 2455087.62 + E x 2.9071 d. 



The RVs thus folded follow a clear double-wave pattern (two 
maxima, two minima). Apart from its obvious double-wave 
shape, the mean radial velocity curve appears slightly anhar- 
monic (the phase difference between the first two periodic sig- 
nals equals 0.34, see Table[3]l. 

3.3. Least-squares deconvolution 

A fit of the mean, observed spectrum of KIC 5988 140 by the 
synthetic one provides quite smooth residuals indicating that 
this is not a double-lined spectroscopic binary. High S/N aver- 
aged profiles compute d by means of the lea st-squares deconvolu- 
tion (LSD) technique dDonati et al. 1119971) confirm this finding. 
Figure [4] illustrates the profiles computed from 1 1 Tautenburg 
spectra (solid) which are shifted in Y-axis for clarity. The phases 
indicated on the plot to the left and to the right of the profiles 
were computed assuming two different periods, Pi and P2, re- 
spectively, and the time of minimum corresponding to the deep- 
est minimum in the Kepler light curve (at BJD~2455085.50, 
cf. Figure [3] bottom part). The LSD-profile computed from the 
mean spectrum which has been used for the detailed analysis in 
Section l3~T1 is shown by a dashed line for comparison. Obviously, 
there is no indication of a secondary component in the spec- 
tra. The first profile from the top appears to be highly asym- 
metric and has a shape typical of that caused by the Rossiter- 



S 




Barycentric Julian Date (2455000.+) 

85 86 87 




Phase (wrt. 2.91 -days period) 



Fig. 3. Top: RVs of KIC 5988140 folded on the period P 2 -2.91 d. 
Hermes RVs are indicated by triangles and TLS ones by open circles. 
Bottom: Portion of the Kepler light curve illustrating the near 2-hr pe- 
riod oscillations. The X-axis shows both BID and phase. In both panels, 
phase zero refers to (B)JD 2455087.62. 
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Fig. 4. LSD-profiles computed based on the individual (solid) and 
mean (dashed) spectra. 



McLaughlin effect (iMcLauglin ll924tlRossiter 11924 . However, 
the computed phases indicate that the corresponding spectrum 
has been taken at the phase near a time of secondary maximum 
light in the Kepler light curve and not during the eclipse phase. 
This fact allows to exclude eclipses as a possible explanation for 
the observed anomalous shape of spectral lines at this phase, thus 
the nature of the anomaly remains unclear. On top of that, there 
are also obvious line profile variations occurring on a short time 
scale which might be caused by stellar oscillations (cf. Sect. [4j. 

4. Kepler photometry 

Similar to the analysis of the spectroscopic data, a search for pe- 
riodicities has been performed using the Period04 programme. 
Figure |5]represents a 2-weeks long string of the Q6 Kepler light 
curve (top) together with the amplitude spectrum based on the 
entire Q6-quarter data string (bottom). The highest peak appears 
at the frequency of 0.687993 d~' (Pi = 1.45350 d), followed by 
the frequency of 0.343984 d _1 (P 2 = 2.90711 d). This result is 
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Table 4. Over 25 frequencies detected in the photometry of 
KIC 5988140. 



6 8 10 

Frequency (d ) 



Fig. 5. Part of the Q6-quarter of the Kepler light curve of 
KIC 5988140 (top) and amplitude spectrum of the entire Q6-quarter 
data string (bottom). 



entirely consistent with the results from the Fourier analysis of 
the RVs (cf. TableO. Note that the phase difference between the 
first two periodic signals is 0.52, therefore not exactly matching 
the one from the radial velocity curve. The ratio F1/F2 equals 
the factor 2, which could in principle be explained as tidal mod- 
ulation in a close binary system (e.g. in ellipsoidal systems), but 
this explanation is not corroborated by the facts that the same 
factor appears in the frequency-analysis of the RVs nor that the 
total RV amplitude is very low. Besides the peaks corresponding 
to the periods Pi and P2, the frequency analysis of the complete 
data set revealed the existence of at least twenty-five more fre- 
quencies, with a significance above ten times the noise level in 
Fourier space. While most of these are located < 0. 1 d _1 indicat- 
ing either long-term changes or remaining instrumental effects 
(see cautionary note from Sect. [2]), some are obviously harmon- 
ics of P2: apart from Pi=P2/2, we also find P2/3, P2/4, and P2/5. 
However, more significant is the detection of nine frequencies 
located beyond 10 d : they all appear to lie in the typical r5 Set 
range (see Table As before, the errors on the frequencies 
and amplitudes were computed based on Monte Carlo simula- 
tions, except for F\ which was kept fixed with respect to F2 (cf. 
Period04). The fraction of the variance removed after prewhiten- 
ing with 27 frequencies (i.e. after removal of the 27 first frequen- 
cies), (1 - R), equals 99%. Figure [3] (bottom) shows a randomly 
selected light curve portion illustrating the short-period varia- 
tions over one and a half full cycle with respect to P2. Besides 
the fact that the light curve appears to be lagging in phase behind 
the RV-curve folded on P2 by at least 0.2 period (i.e. maximum 
to maximum), a very clear beating of the most significant fre- 
quency located at 14.32981 d -1 (corresponding to a pulsation 
period of lh40m) can be seen, caused by the presence of six 
more frequencies located between 14 and 15 d . Among these, 
a few frequency spacings are repeatedly occurring. For example, 
the frequencies F3, Fj and F23 are separated by almost exactly 
0.040 d -1 . The couples (F5,Fg) and (F\g,Fs) show a difference of 
0.31 and 0.25 d~ 1 , respectively. Two other <5 Set-type frequencies 
are F A = 12.28328 d" 1 and F 14 = 10.91996 d l (with respective 
period ratios of 0.86 and 0.76). The latter is an interesting pe- 
riod ratio, close to the ratio 1H/F expected for radial modes in 
such a star. The occurrence of Fn = 23.89006 d -1 , on the other 



ID 


Frequency 


Amplit. 


S/N 


Note 




(± error) 


(± 0.001) 








d 1 


ppt 


ppt 




F t 


0.687993 (fixed) 


1.261 


411 


= 2F 2 


F 2 


0.3439842 (3E-07) 


1.123 


335 




F 3 


14.32981 (2E-06) 


0.910 


1207 




F4 


12.28328 (4E-06) 


0.192 


256 




F s 


14.64803 (5E-06) 


0.146 


230 




Fi 


14.28878 (9E-06) 


0.094 


123 




Fg 


14.33072 (2E-05) 


0.084 


112 




F10 


1.031966 (1E-05) 


0.065 


33 


= 3F 2 


1 1 5 


14.12919 (1E-05) 


0.055 


70 




Fid 

* 14 


10.91996 (2E-05) 


0.051 


129 




F \ 7 


23.89006 (2E-05) 


0.045 


49 




* 1 8 


1.375905 (2E-05) 


0.044 


38 


= AFn 


^ \9 


14.90637 (2E-05) 


0.043 


77 




F 9 q 


14.24855 (2E-05) 


0.033 


43 






1.719992 (3E-05) 


0.028 


33 


= 5F 2 


n 


0.057054 (1E-05) 


0.156 


45 




n 


0.053534 (1E-05) 


0.092 


26 






0.060134 (2E-05) 


0.094 


27 






0.090642 (1E-05) 


0.066 


19 




1 J 


0.075828 (1E-05) 


0.065 


19 






0.083162 (2E-05) 


0.049 


14 




F 


0.004253 (2E-05) 


0.044 


13 






0.033734 (2E-05) 


0.036 


20 




1 


0.062921 (2E-05) 


0.048 


14 




F u 


0.050308 (2E-05) 


0.042 


12 




F'25 


0.092842 (2E-05) 


0.036 


10 






0.070695 (3E-05) 


0.035 


10 






0.110442 (3E-05) 


0.029 


8 





Frequencies not considered as meaningful, see Sect. [2] 



hand, is questionable, since it lies very close to the Nyquist fre- 
quency of the data set. We further notice the presence of various 
long-term periodicities (frequencies < 0.01 d -1 ) which appear 
with a low significance, though none as small as the frequency 
of 0.001 d _I found below the detection limit in the RV data (cf. 
Sect. I3.2l i. We also recall that such very low frequencies might 
be affected by the data (pre)processing algorithms. 

Next, we tested the stability of the properties of the <5 Set- 
type pulsations by searching for a best-fit model allowing for 
a periodic variation of the phase shifts of these frequencies in 
the multi-parameter solution. For this, we used the model called 
'Periodic Time Shift' in Period04 (cf. PTS model under expert 
mode). The computations show that no better fit can be found 
by including a possible (long-term) variation of the phase shifts 
which might possibly be due to a modulation by either one of the 
low frequencies Pi or P2. Our conclusion is that the properties 
of the 6 Set-type frequencies are stable. 

We also verified the stability of the shape of the "basic" 
light curve (i.e. cleaned for the near 2hr-period variations and 
the long-term trends) based on the main frequency and the har- 
monics by prewhitening with all other frequencies detected dur- 
ing the Fourier analysis, and by plotting the folded residual light 
curve using colours for different subsets. Figure [6] shows the re- 
sulting phase diagram when subdividing the data into bins of 
tens of days. We see no obvious temporal change in the shape of 
this curve. As an independent test, we divided the "basic" Kepler 
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Fig. 6. "Basic" light curve of KIC 5988140 based on all the data and 
folded on the period P 2 (after prewhitening for 23 frequencies). 



light curve into different cycles and fitted each with a 4th-order 
trigonometric polynomial according to 



m = A + ^] Aj ■ sin {in ■ jF 2 ■ t + cf>j) , 



(1) 



7=1 



where m is the magnitude, Aj is the amplitude, jF 2 stands for the 
frequency F2 and its harmonics, t is the time of the observation, 
<f>j is the phase and index j runs from 1 to 4. Then, we char- 
acteri sed each cycle by its Fourier coefficients dSimon & Teavs I 
1 19821) . In Fig. [7] we represent the amplitude, A\, and the rela- 
tive variations of the parameters /?2i = A2/A1, = A3/A1 as 
a function of time. The amplitude A\ shows a scatter of about 
0.08 ppt which is ten times larger than the average fit error 
(~0.008 ppt). While the relative scatter of R 2 \ is small, the rela- 
tive scatter of Rt,\ is more important, but considering that the am- 
plitude A3 is much smaller, than A\, the corresponding variation 
remains insignificant and we find no clear trend. This demon- 
strates that the "basic" light curve is stable over a period of al- 
most two years (T = 682 days). 

Considering that the time resolution in Long Cadence mode 
is of the order of 30 min, we must caution the determination of 
the amplitudes of the 6 Set-type frequencies in Table[4] Whereas 
the amplitudes in the low-frequency regime may be correctly 
evaluated, the ones in the high-frequency regime might actually 
be slightly underestimated. 

Making use of Figures [3] and [6j we can also verify the 
phase relationship between the "basic" light and RV curves: both 
curves are not exactly in anti-phase with respect to each other: 
the moments of minimum RV precede the moments of maximum 
light by about 0. 1 period. 

After prewhitening for all the frequencies listed in Table |U 
the final residuals (see Fig. [8]) do not look perfectly flat, as one 
would expect, but reflect small instrumental effects and reveal 
discontinuities which were not accounted for during the homo- 
geneization (pre)processing. 



5. Discussion 



It seems that lUvtterhoeven et al~l d201 lbl) based their binary 
classification on the morphology and the high regularity of the 
light curve (see Figure [5]), where the deepest minima, occur- 
ring every ~2.91 days, could be attributed to primary eclipses. 
Our Fourier analysis of the Kepler data reveals the existence 
of multiple frequencies, with the highest amplitude peaks oc- 
curring at the frequencies of 0.687993 d '(Pi = 1.45350 d) 
and 0.343984 d _1 (P 2 = 2.90711 d), with an exact integer ra- 
tio of 2, the first harmonic term being the most significant. 
A Fourier analysis of the ground-based RV data indicates the 
same frequencies, moreover in the same order of importance, 
hence showing the double-wave shape. The fact that a double- 
wave pattern, almost but not exactly in anti-phase with the light 




5200 5400 
BJD (2450000+) 



Fig. 7. Temporal evolution of the (relative) amplitudes A\ and , ^31 
of the "basic" light curve of KIC 5988140. 





Fig. 8. Residual light curve of KIC 5988 140 based on all the data (after 
prewhitening for 28 frequencies). 



curve (Sect. @}, is also detected in the radial velocity phase dia- 
gram excludes binarity as a possible cause of the observed low- 
frequency variability. Moreover, the total RV amplitude of al- 
most 8 km s is too small for a short-period binary system with 
a primary of spectral type A, unless the inclination angle would 
be very low (in that case, we would not detect any feature resem- 
bling an eclipse). 

We next consider the model of rotational modulation gener- 
ated by an asymmetric intensity distribution on the stellar sur- 
face. The latter could be (a) due to the presence of stellar spots 
on the surface, or (b) caused by a (thin) convection zone located 
on or near the surface. The RVs folded on the period of ~2.91 d 
(see Figure|3] top) show a double wave which could be explained 
by a surface structure that includes two spots diametrically op- 
posed to one another. The same is also true for the double-wave 
pattern detected in the Kepler light curve. That a harmonic of 
the rotation period of the star has the highest amplitude in the 
RV data, was previously found i n some spotted stars, e .g., in the 
He- weak silicon star HR7224 dLehmann et al. 1120061) . On the 
other hand, both the spectral type (late A) and the long-term sta- 
bility of the "basic" light curve over a period of about 2 yrs argue 
against variability due to stellar spots for this object. 

Usually, we don't expect magnetic activity (causing spots 
on the stellar surface) nor convection to occur in the atmo- 
spheres of normal A-type stars. However, convection does oc- 
cur in some late A-type stars, and particularly in the photo- 
sphere where it competes with radiation transporting the flux. 
The case of Altair is illustrative: this fast rotating A-type star 
has a broad equatorial band at T e +T = 6900 K which gives 
rise to str ong convection and chromospheric activity like C II 
emissi on dPeterson et al.ll2006|). Us ing an XMM-Newton obser- 
vation, iRobrade & Schmittl d2009h confirmed the presence of 
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(coronal) X-ray emission and weak magnetic activity. Altair in- 
deed belongs to the transition region where stars can develop 
envelopes in-between fully radiative and fully convective ones. 
Based on observations co llected by the Wide Field Infrared 
Explorer (WIRE) satellite, iBuzasi et al~l d2005l) showed that it 
is also a very low-amplitude 6 Set star (with Am < 1 ppt). 
Interestingly, in addition to the frequencies detected in the 15- 
29 d _1 range, the existence of two frequencies in the low range, 
namely at 3.53 and 2.57 d -1 (cf. thei r Table 1) is demonstr ated. 
Using long-baseline interferometry, iPeterson et al. I d2006h de- 
rived the rotational frequency of 2.71 d" 1 , which i s close to the 
lowe st frequency detecte d by IBuzasi et al. I d2005l). I n essence, 
both lOhishi et al~l (12004 and IPeterson et al. I (l2006h conclude 
that the surface of Altair displays an extremely asymmetric in- 
tensity distribution and that the asymmetry is consistent with 
that expected from the known high rotation and oblateness. It 
thus appears that the rotationally distorted stellar surface of 
Altair induces low amplitude modulations which are detectable 
in the high-quality light curves from space. We can ask our- 
selves whether KIC 5988140 might be similar to Altair or to 
the weakly active, X-ray emitting A5/F0-type star HR 8799 
dRobrade & Schmittll2010l) . 

Adopting the radius of 3.57 R Q (/?») listed in the KIC, we 
estimate a rotation period of ~3.7 d (reciprocal of 0.27 d~') and 
of ~2.5 d (reciprocal of 0.40 d~') for an inclination angle of 
90° and 45°, respectively. The detected period of 2.9071 1 d lies 
in-between these two derived periods. We point out that there 
is also a triplet of frequencies (^3,4,5) which has a frequency 
spacing of the order of 0.3 d~' that could be explained as 
rotational splitting of a single, independent frequency. Using 
V e quat = 50.61 -R*/P2 = 62.3 km s _1 and the measured v sin i we 
obtain an inclination angle of ~50°, which is a most reasonable 
value. We next attempted to find a model including two hot 
spots symmetrically located on opposite sides of the stellar 
equator in order to explain the patterns in both the light and 
the RV curves. For this, we used a simple model in which the 
stellar surface has intensity /o, including a standard law for the 
limb darkening, the spots have a size FWHM, (corresponding 
to the Full Width at Half Maximum or FWHM) and a Gaussian 
intensity distribution of amplitude A,. The following formula 
was adopted: 

///o = (1 + spot\ + spoti) • limb 
spoU = Ai ■ exp [- \n(2)(ri/FWHMi) 2 ] 

and 

limb = 0.4 + 0.6 ■ y/l-z 2 , 

where r, is the distance on the stellar surface from the centre of 
spot ; and z is the distance from the centre of the visible disk. For 
the intrinsic absorption line profile, we used a Gaussian having a 
width corresponding to £ = 3.16 kms -1 . The resulting line pro- 
files are obtained by integrating over the visible surface where 
each point has a velocity corresponding to vsin/ = 52 kms -1 . 
Deduced RVs represent the first moments of the line profiles. 

We chose spot sizes and amplitudes such as to reproduce the 
observed total amplitude of the RV curve. Figure [9] illustrates 
two cases. In the first case (model A), we chose a maximum size 
of both spots equal to the stellar radius. An amplitude A i = 0.78 
was needed to reproduce the observed RV amplitude of 8 km s _I . 
In the second case (model B), we chose a size FWHM = R*/30 
and needed an amplitude A i = 23 (A 2 was chosen a bit smaller in 
each case to reproduce the observed asymmetry). Model B was 
(only) selected as the opposite limit case, since such a small and 



Table 5. Results of the analysis of four individual spectra of 
KIC 5988140 taken at four different phases of the period P 2 . 
T e ff is in K, vsin; and £ in kms 1 , \ogg, [M/H] and the indi- 
vidual abundances [A/H] are in dex. The uncertainties are listed 
in terms of last digits in parentheses. 



Param 


Spectr. 1 


Spectr. 2 


Spectr. 3 


Spectr. 4 


Cm-. 

oiin 






Fundamental parameters 






T c ff 


7680(50) 


7590(50) 


7630(50) 


7640(50) 




logg 


3.58(15) 


3.47(15) 


3.51(15) 


3.57(15) 




v sin/ 


51.0(1.5) 


51.5(1.5) 


50.5(2.0) 


52.0(2.0) 






3.12(20) 


2.95(25) 


3.30(20) 


3.35(25) 




[M/H] 


-0.28(05) 


-0.31(05) 


-0.30(05) 


-0.30(05) 








Individual abundances 






[Fe/H] 


-4.89(05) 


^1.94(05) 


^1.99(05) 


-4.99(05) 


-4.59 


[Ca/H] 


-6.03(10) 


-6.03(10) 


-6.08(10) 


-6.13(10) 


-5.73 


[Ti/H] 


-7.44(10) 


-7.44(10) 


-7.44(10) 


-7.44(10) 


-7.14 


[Cr/H] 


-6.60(15) 


-6.65(15) 


-6.65(15) 


-6.65(15) 


-6.40 


[Mg/H] 


-4.61(15) 


-4.61(15) 


-4.61(15) 


-4.66(15) 


-4.51 


[Ni/H] 


-5.96(15) 


-6.01(15) 


-6.01(15) 


-6.01(15) 


-5.81 


[Sc/H] 


-9.09(20) 


-9.19(20) 


-9.14(20) 


-9.09(20) 


-8.99 


[Mn/H] 


-7.05(20) 


-7.10(20) 


-7.10(20) 


-7.10(20) 


-6.65 


[V/H] 


-8.34(25) 


-8.39(25) 


-8.39(25) 


-8.34(25) 


-8.04 


[C/H] 


-3.75(25) 


-3.85(25) 


-3.80(25) 


-3.75(25) 


-3.65 


[Y/H] 


-9.93(30) 


-10.03(30) 


-10.03(30) 


-9.93(30) 


-9.83 


[Si/H] 


^1.68(30) 


-4.78(30) 


-4.68(30) 


-4.73(30) 


-4.53 


[Co/H] 


-7.47(35) 


-7.42(35) 


-7.37(35) 


-7.42(35) 


-7.12 


[Sr/H] 


-8.92(40) 


-8.97(40) 


-9.12(40) 


-9.17(40) 


-9.12 


[Ba/H] 


-9.27(50) 


-9.47(50) 


-9.47(50) 


-9.37(50) 


-9.87 


[Eu/H] 


-11.52(50) 


-11.42(50) 


-11.47(50) 


-11.57(50) 


-11.52 



bright spot would produce bumps moving across the line pro- 
files (which were not detected). Figure [10] shows the resulting 
RV and light variations. Both models predict intensity changes 
with a total amplitude of about 20%, i.e. a factor of 40 larger than 
the observed one (max. 5 ppt peak-to-peak). Varying the shape 
of the spots (exponential instead of Gaussian), their position in 
latitude, their separation in longitude, or the inclination of the 
rotation axis did not yield a better agreement. Thus, the small 
value of the observed light-to-velocity amplitude ratio cannot be 
explained by a (simple) spotted surface model. The same argu- 
ment also invalidates the assumption that KIC 5988 140 might 
have a partially convective atmosphere (along with the fact that 
it has a low probability of being a fast rotator with v sin i = 
52 km s^'assuming that P2 is the rotation period). 

Spots will show as regions of different chemical composition 
on the stellar surface. To verify this, we performed separate anal- 
yses of four Hermes spectra taken at four different phases of the 
period P2 (e.g. the extrema of the light curve) to look for possi- 
ble changes in the derived abundances and/or stellar parameters. 
These results are summarized in Table [5] The largest deviations 
in T e ff, log g, and £ are respectively 80 K, 0. 1 dex and 0.3 km s _1 , 
which is comparable to the uncertainties listed in Table [2] Both 
the metallicity and v sin i are stable. We may therefore conclude 
that there is no obvious difference in the obtained fundamen- 
tal parameters nor in the individual abundances with respect to 
phase. 

Third, nine of the photometrically detected frequencies 
(-^3,4,5' Fj, Fg i io ) F13, F19, Fn, see Table |4]i lie in the typical 
6 Set range. Observed regular spacings are of the order of 0.3 and 
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Fig. 9. Spot models A (left) and B (right) shown for a certain rota- 
tion phase and drawn on separate greyscales to cover the full brightness 
range of both cases. 
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Fig. 10. RV and relative intensity variations for spot models A (solid) 
and B (dashed). 



0.04 cT 1 . The question thus remains if the frequencies F\ and F2 
(and the harmonics) can possibly be related to the y Dor-like 
oscillations, which would confirm its class ification as a 6 Sct- 
y Dor hybrid star bv lCatanzaro et al."1 d201 ll) . However, and most 
importantly, the same reasoning as presented in the model in- 
volving rotational modulation is applicable here: if the RV vari- 
ations were caused by pulsations, we would expect to detect light 
variations with a total amplitude in scale with the one of the RVs 
(assuming that low-degree modes only are excited). The light-to- 
velocity amplitude ratio for slowly pulsating B stars lies around 
4-5 mmag/ kms 1 (based on 13 con firmed cases, cf. top panel of 
Fig. 19 in iDe Cat & Aerts 1120021) . Assumin g that this ratio is 
also valid for y Dor stars (e.g. from Table 2 in lAerts et al. 1120041 
, we have a mean amplitude ratio of 15 mmag/kms -1 from two 
y Dor stars), we would expect a total amplitude of about 16-20 
mmag in the "basic" Kepler light curve, i.e. 15-20% (similar to 
the previous scenario). We are (again) several orders below this 
expected value with the observed total amplitude of the "basic" 
light curve (max. 5 ppt). Indeed, the problem of the unbalanced 
light-to-velocity amplitude ratio remains with this explanation 
too. Other, less critical, counter-indications are: (1) y Dor stars 
are typically multi-periodic pulsators; (2) the small integer ratio 
of 2 between the two most dominant frequencies is unexpected 
in the case of (non-radial) pulsators which are not a member of a 
binary system, (3) the existence of more harmonic frequencies, 
(4) this star is perhaps too evolved to be a y Dor-type pulsator (cf. 
its position in the H-R diagram), and (5) the previously derived 
phase relationship (Sect.|4| is unusual in the case of pulsation, as 
the moments of minimum RV (i.e. maximum surface layer ex- 



pansion) generally follow the moments of lig ht maximum (e.g. 
in 6 Scuti stars, see iPreston & Landolt1ll99S 



6. Conclusion 

KIC 5988 140 (HD 188774) was assigned a spectral type of 
A7.5 IV-III which, along with its atmospheric parameters and 
abundances (see Tables Q] and |2|, puts it in the middle of the 
6 Set and just at t he blue edge of the y Dor instability strips 
dFekel et al. 112003b . The presence of the multiple frequencies in 
the range 10-15 d _1 detected from the Kepler photometry is firm 
evidence that the star is a 5 Set pulsating star. The apparently 
slightly more dominant low frequencies detected in both the RV 
data and in the Kepler photometry, identified as one major fre- 
quency with harmonics, may be caused by (at least) three very 
different physical processes: a) binarity, b) rotation in combina- 
tion with an inhomogeneous surface intensity/chemical distribu- 
tion or c) y Dor-type pulsation in the gravity-mode regime. 

In our case, we are left with an open question regarding 
the phenomenon which causes the detected low frequencies in 
KIC 5988 140 as none of the abovementioned scenarios gives an 
entirely satisfactory explanation for the observed variability pat- 
terns discovered in this late A-type star. It seems indeed difficult 
to reconcile the Kepler light curve with the ground-based RV 
curve. We discarded binarity on the basis of the double-wave 
pattern detected in the RV curve. This is a strong argument. The 
light-to-velocity amplitude ratio is an obstacle for the other two 
scenarios: that of rotational modulation as well as that of pul- 
sations. It might be that our simple modelling of stellar spots is 
lacking some reality, but we believe that it, at least qualitatively, 
shows the order of the expected effects. 

This reminds us of another still unexplained behaviour ob- 
served by Kepler, the fact that the Kepler light curves of many 
A- an d F-type stars show variability at low frequencies dBalona I 
1201 ll) . The author analysed the distributions of v sin i derived on 
the basis of the (dominant) low frequencies and those of the ob- 
served rotational frequencies and showed that both agree rather 
well assuming that the value of the dominant low frequency is 
the rotational period. In the specific case of the 6 Set stars ob- 
served by Kepler, assuming that half of t he value of the d ominant 
low frequency is the rotational period, iBalonal d201 ll) found a 
good agreement between bot h (observed and derived) v sin i dis- 
tributions. On the other hand. lGrigahcene et al. I (120101) proposed 
that most Kepler 8 Set stars are of the hybrid type where the low 
frequencies are due to the y Dor phenomenon. An explanation 
for this behaviour is currently needed. 

KIC 5988 140 is most likely not unique amongst the Kepler 
A-type stars. However, it is one of the first cases where a high- 
quality data set in RV has been collected and confronted to the 
photometry. Because of the difficulty in finding a plausible ex- 
planation for the cause of the observed variations, we might con- 
sider the more complex model of a triple system to explain the 
double wave detected in the RV curve. We can exclude a hierar- 
chical system consisting of a single star and a binary based on 
the stability criterion because of the 1 :2 ratio of the long-term 
frequencies. On the other hand, a single star revolved by two 
satellites could be stable. From the periods and the RV ampli- 
tudes, using a typical mass of 2-3 M for the A-type star, we es- 
timate masses in the range of brown dwarfs to M-dwarfs for both 
objects, namely in the range 20-220 Jupiter masses depending 
on the orbital inclination angle. However, these satellites would 
have very small semi-major axes of about 8 and 12 R Q and we 
cannot tell anything yet about the stability of such close orbits in 
a 1:2 resonance. From stability considerations, the low-mass ob- 
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jects (i.e. brown dwarfs) seem to be more likely. Whether such a 
scenario might also explain the observed photometric variations 
by reflection (including ellipsoidal variability as well as beam- 
ing effects), is presently uncertain and will require an in-depth 
analysis. We present this idea merely as a possible working hy- 
pothesis for our future work on this intriguing Kepler star. 

One should pay attention to other cases in the Kepler field 
which might show the same difficulty in interpretation. This re- 
minds us that, without complementary (even ground-based) in- 
formation in the form of spectroscopy and/or multi-colour pho- 
tometry, any interpretation based on single-passband light curves 
alone may be biased. A valuable method allowing to distin- 
guish between binarity rotational modulation or pulsations of 
type y Dor would be to collect a mplitude ratio's in multiple pass- 
bands (e.g. iHenrv et al. 1120071) . To perform such a task at the 
required scale and accuracy level will however most certainly 
necessitate future space resources. 

As illustrated by this study, the acquisition of high-quality 
ground-based spectra can provide relevant, new information for 
a more thorough understanding of the hyper-quality white-light 
stellar photometry provided by the space missions. 
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Appendix A: Table A.1 to be published 
electronically 



Table A.l. Journal of recent spectroscopic observations. Listed are 
the sequence number, the exposure time, the barycentric Julian Date, 
the measured RV, and the accuracy in RV. 



Nr 


xp 


Date 


BJD 


RV 


Error RV 




_tsec)_ 




(2455000+) (kms- 1 ) (kms-') 








Hermes 






1 


1200 


2010/06/05 


3^3 704S428 


-18.942 


0.020 


2 


1800 


2010/07/14 


392.5180032 


-19.627 


0.032 


3 


1800 


2010/07/14 


392.6444271 


-19.964 


0.021 


4 


2000 


2010/07/16 


394.5169717 


-18.502 


0.020 


5 


1800 


2010/07/16 


394.6163039 


-17.147 


0.025 


o 


9700 
Z 1 UU 


901 0/07/1 8 
ZU 1U/U//10 


396.4641905 


-13.775 


0.021 


7 


9AOO 
Z<+UU 


901 0/0771 8 
ZUtU/U//to 


396.6385025 


-16.191 


0.029 


o 


9700 
Z 1 UU 


901 0/07/90 
ZU 1U/U //zu 


398.4383502 


-20.191 


0.021 


V 


9ztOO 
ZH-UU 


901 0/07/90 
ZU 1U/U //zu 


398.6246792 


-19.895 


0.029 


i u 


1 900 
1 zuu 


901 O/OS/O^ 
ZUlU/Uo/UJ> 


412.4096963 


-15.793 


0.019 


1 1 


^00 
JUU 


901 O/OS/O^ 
ZUlU/Uo/UJ 


414.4434733 


-19.671 


0.013 


1 9 

1 z 


^OO 
jUU 


901 O/OS/OS 
ZUlU/Uo/UJ 


414.6630946 


-20.497 


0.035 




^OO 
jUU 


901 O/OS/07 
ZU 1U/UO/U / 


416.4173718 


-14.901 


0.018 


1 A 


^00 
JUU 


901 0/0S/07 
ZU 1U/UO/U / 


416.7144405 


-13.203 


0.031 


1 K 

1 J 


^00 
JUU 


901 O/OS/OQ 

zuiu/uo/uy 


418.3981617 


-18.337 


0.034 


16 


300 


2010/08/09 


418.7080667 


-20.417 


0.026 


17 


300 


2010/08/11 


420.4033010 


-20.640 


0.032 


18 


300 


2010/08/12 


421.4271506 


-18.821 


0.033 


19 


600 


2010/08/12 


421.5680935 


-19.948 


0.029 


20 


300 


2010/08/14 


423.4211612 


-20.538 


0.025 


21 


480 


2010/08/14 


423.5536145 


-19.015 


0.014 


22 


600 


2010/08/15 


424.4009403 


-19.750 


0.020 


23 


600 


2010/08/16 


425.3884691 


-13.305 


0.030 


24 


600 


2010/08/16 


425.5342909 


-14.257 


0.023 


25 


550 


2010/08/17 


426.3884962 


-20.266 


0.028 


26 


550 


2010/08/17 


426.6004571 


-17.948 


0.015 


27 


600 


2010/08/18 


427.4912895 


-20.324 


0.028 


28 


400 


2010/08/18 


427.5044292 


-20.821 


0.029 


29 


1800 


2010/08/20 


429.5096819 


-17.128 


0.031 








TLS 






30 


900 


2011/07/12 


754.5635120 


-19.917 


0.022 


31 


1800 


2011/07/14 


757.4733430 


-20.035 


0.027 


32 


1800 


2011/08/08 


782.3464050 


-18.128 


0.033 


33 


1800 


2011/09/10 


815.4850340 


-18.312 


0.028 


34 


1800 


2011/09/11 


816.4960600 


-15.304 


0.030 


35 


1024 


2011/09/12 


817.2843500 


-18.059 


0.032 


36 


1800 


2011/09/12 


817.4911610 


-14.930 


0.045 


37 


1800 


2011/09/13 


818.3643670 


-18.410 


0.028 


38 


1800 


2011/09/14 


819.2709620 


-14.222 


0.020 


39 


1800 


2011/09/14 


819.2925580 


-14.563 


0.025 


40 


1800 


2011/09/19 


824.4577250 


-19.448 


0.033 
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